Finite temperature excitations of a trapped Bose-Fermi mixture 
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We present a detailed study of the low-lying collective excitations of a spherically trapped Bose- 
Fermi mixture at finite temperature in the collisionless regime. The excitation frequencies of the 
condensate are calculated self-consistently using the static Hartree-Fock-Bogoliubov theory within 
the Popov approximation. The frequency shifts and damping rates due to the coupled dynamics of 
the condensate, noncondensate, and degenerate Fermi gas are also taken into account by means of 
the random phase approximation and linear response theory. In our treatment, the dipole excitation 
remains close to the bare trapping frequency for all temperatures considered, and thus is consistent 
with the generalized Kohn theorem. We discuss in some detail the behavior of monopole and 
quadrupole excitations as a function of the Bose-Fermi coupling. At nonzero temperatures we find 
that, as the mixture moves towards spatial separation with increasing Bose-Fermi coupling, the 
damping rate of the monopole (quadrupole) excitation increases (decreases). This provides us a 
useful signature to identify the phase transition of spatial separation. 

PACS numbers: PACS numbers:03.75.Kk, 03.75.Ss, 67.60.-g, 67.40.Db 
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I. INTRODUCTION 

The impressive experimental achievement of Bosc- 
Einstein condensation (BEC) in the bosonic systems 
87 Rb 0, 23 Na 0, and 7 Li has initiated and stimulated 
a whole new field of research in the physics of quantum 
atomic gases. Recently, several groups have extended 
these experiments to the case of trapped Bose-Fermi mix- 
tures, in order to employ the "sympathetic cooling" to 
reach the regime of quantum degeneracy for the Fermi 
gas. As a first step, stable Bose-Einstein condensates 
immersed in a degenerate Fermi gas have been realized 
with 7 Li in 6 Li 4), 23 Na in 6 Li 5], and very recently 



with 87 Rb in 40 K [||. In the possible next step, inves- 
tigations of the thermodynamics, collective many-body 
effects and other properties could be available soon in 
these systems. Especially interesting is the behavior of 
low-energy collective excitations, since the high accuracy 
of frequency measurements and the sensitivity of collec- 
tive phenomena to the interatomic interaction make them 
good candidates to unravel the dynamical correlation of 
the many-body system. 

On the theoretical side, several analyses have been pre- 
sented for low-lying collective excitations of a trapped 
Bose-Fermi mixture. Collective modes in the collisionless 
limit, where the collision rate is small compared with the 
frequencies of particle motion in traps, have been con- 
sidered by the sum-rule approach , by the scaling the- 
ory , or in the random-phase approximation 0, ll(| . 
In the collision-dominated regime collective oscillations 
have been discussed by Tosi et al. in Refs. 0,0, and 
by the authors in Ref. 0. These investigations have 
mainly concentrated at zero temperature using the stan- 
dard two-fluid model for the condensate and the degener- 
ate Fermi gas. However, the realistic experiment is most 



likely carried out at relatively higher temperatures, where 
the condensate oscillates in the presence of a considerably 
large fraction of above-condensate atoms. It thus seems 
timely to develop an extension of these theories to the 
finite temperature. 

In the present paper we investigate the low-lying collec- 
tive excitations of a spherically trapped Bose-Fermi mix- 
ture at finite temperature in the collisionless regime. We 
confine ourselves to the collective modes of the conden- 
sate, i.e., the density oscillations of the condensate. We 
first calculate the mode frequencies by using the simplest 
temperature-dependent mean-field theory — the Popov 
version of the Hartree-Fock-Bogoliubov (HFB) theory 
- that has been generalized by us to a tra ppe d Bose- 
Fermi mixture to study its thermodynamics [13j. For a 
purely Bose gas, it is well known that the HFB-Popov 
theory includes only the static mean-field effects of the 
noncondensate atoms [l4| and thus predicts the correct 
mode frequency only at temperatures T < To = 0.65T C 
[Tlj . where T c is the critical temperature of BEC. Above 
To the noncondensate component becomes considerably 
large and its dynamics should be treated on an equal 
footing with that of the condensate [H HI 111 ITl l20| . 
In our case of Bose-Fermi mixtures, the situation is more 
crucial. Due to the large number of fermions, the coupled 
dynamics of the condensate, noncondensate, and degen- 
erate Fermi gas has to be taken into account even at 
zero temperature. In this paper, we shall treat it per- 
turbatively in the spirit of the random phase approxima- 
tion (RPA) and linear response theory. We derive the 
explicit expression for the frequency shift and damping 
rate arising from the coupled dynamics, which in the ab- 
sence of the Bose-Fermi interaction coincides with the 
finding of Ref. [0j- Based on this expression and the 
static HFB-Popov theory, we present a detailed numer- 



2 



ical study of the monopole and quadrupole condensate 
oscillations against the Bose-Fermi coupling. The dipole 
excitation is also studied and found to be consistent with 
the generalized Kohn theorem. 

The paper is organized as follows. In the next sec- 
tion we derive the theory used in this paper. In Sec. 
Ill we apply this theory to mixtures of 41 K— 40 K and 
87 Rb— 40 K, and calculate the dispersion relation of the 
monopole and quadrupole excitations as a function of 
the Bose-Fermi coupling. The behavior of monopole and 
quadrupole modes against temperature is also discussed 
in detail. Finally, section IV is devoted to conclusions. 



II. FORMULATION 

In this section we first generalize a time-dependent 
mean-field scheme developed by Giorgini ^(| to Bose- 
Fermi mixtures, and derive the equation for the small- 
amplitude oscillations of the condensate. Since the for- 
malism of this time-dependent mean-field approximation 
for an inhomogeneous interacting Bose gas has already 
been presented in detail in Ref. |l6| , here we shall merely 
concentrate on the key points, and indicate the necessary 
modification in the presence of the fermionic component. 
By means of the RPA and linear response theory, we fur- 
ther consider the fluctuations of the noncondensate and 
of the degenerate Fermi gas induced by the condensate 
oscillations. The back action of these fluctuations on 
the condensate motion is then calculated perturbatively 
to second order in the interaction coupling constant to 
obtain the explicit expression for frequency shifts and 
damping rates. 

Our starting point is the trapped binary Bose-Fermi 
mixture that is portrayed as a thermodynamic equilib- 
rium system under the grand canonical ensemble whose 
thermodynamic variables are N b and Nf, respectively, 
the total number of trapped bosonic and fermionic atoms, 
T, the absolute temperature, and [i b and /i/, the chemi- 
cal potentials. In terms of the creation and annihilation 
bosonic (fermionic) field operators ip + (r,t) and ip( r it) 
(</> + (r, t) and 4>(r, t)), the density Hamiltonian of the sys- 
tem takes the form (in units of Ti = 1, and all field oper- 
ators depend on r and t) 
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(1) 



Here we consider a spherically symmetric system, with 
static external potentials V^l p (r) — m b jcu^r 2 /2, where 
mbj are the atomic masses, and u> b j are the trap fre- 
quencies. The interaction between bosons and between 
bosons and fermions are described by the contact po- 
tentials and are parameterized by the coupling constants 



cjbb = 4:irh~ abb /rrib an d 9b f — 2itTi a b f/m r , to the lowest 
order in the s-wave scattering length a bb and a b f, with 
m r = m b mf/(m b +rrif) being the reduced mass. 



A. time-dependent mean-field approximation 

According to the usual treatment for Bose system with 
broken gauge symmetry, we shall apply the decomposi- 
tion: ip(r,t) — 3>(r, t) +-0(r,i), where $(r,i) = (ip(r,t)) 
represents a time-dependent condensate wave function 
and allows us to describe situations in which the system is 
displaced from equilibrium and the condensate is oscillat- 
ing in time. With this respect the average (...) is intended 
to be a non-equilibrium average, while time-independent 
equilibrium averages will be indicated in this paper with 
the symbol (...) . The field operator ip(r,t) plays the 
role of excitations out of the condensate, and by defini- 
tion satisfies the condition ^f/j(r,t)^) = 0. This ansatz is 
then inserted in the equation of motion for ip(r, t): 
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V 



Ivrib 

+9bbip + ^ip + gbfip<t> + <t>- (2) 

Taking a statistical average over Eq. (|2J) and setting 
the triplet average values (ip + (r, i)V>(r, t)t/>(r, t)\ and 

(^(r,i)<?(> + (r,i)0(r,t)^) to zero H2 thus leads to the fol- 
lowing equation of motion for the condensate wave func- 
tion 



+g bb |$(r, t)r $(r, t) + 2g bb n(r, i)$(r, t) 
+<? 6b m(r, (r, t) + g b fn f {r, i)$(r, i)(3) 

where the densities are defined, respectively, as h(r,t) = 
^+(r,t)^(r,t)V m(r,t) = /^(r, f)V>(r, f)^ , and 

n/(r, t) = (4> + (r, t)<fr(r, t)). Under the stationary con- 
dition, we replace <I>(r,i), h(r, i), m(r, i) and n/(r, t) 
by their equilibrium values $o( r ) = (V ? ( r ))o) ^°( r ) = 



ip + '{r)ip(r)J , m°(r) = ^-0(r)-0(r)y and n°(r) = 

{(f> + (r)(j)(r)) , respectively. This yields the general- 
ized time-independent Gross-Pitaevskii (GP) equation 
for Bose-Fermi mixtures |l3| 



V 

2m b 



V f b 

y trap 



fJ-b + 9bb (n + 2n°) + g b fn°j 



$o = 0, 



(4) 

where no(r) — |$o( r )| is the condensate density. In the 
above equation, we already use the Popov prescription: 
to (r) = 0, which amounts to neglect the effects arising 
from the equilibrium anomalous density |25| . 

We are interested in the small amplitude oscillations of 
the condensate, which is only slightly displaced from its 



3 



stationary value $o( r ) : &( r , t) — < &o(r) + S$(r, t), where 
<5<I>(r, t) is a small fluctuation. This small oscillations can 
consequently induce small fluctuations of the densities 
around their equilibrium values: n(r,t) = h (r)+5n(r,t), 
m(r, t) — 5m(r, t), and n/(r, t) = n° (r) + <5n/(r, t). The 
time-dependent equation for <5$(r,f) is then obtained by 
linearizing the equation of motion Q3J 



£J$(r,t)+5»no(r)*$*(r,t) 

+2 5b6 $ (r)fe(r, t) + .g 6b $ (r)(5m(r, t) 
+# 6/ $o(r)<5n/(r,t), (5) 

where we have introduced the Hermitian operator 



C = - 



V 
2mj 



+ ^trop (r) - + 2g bb n b (r) + g bf n° f (r), (6) 



and n° (r) = rio(r) + n°(r), the total density of bosons. 
In Eq. (jSJ) , the terms containing 5h, 8fh and 5n / account 
for the dynamic coupling between the condensate and 
the fluctuations of the noncondensate component and of 
the degenerate Fermi gas. Assuming that the condensate 
oscillates with frequency u: S$(r,t) — <5$(r)e~ 4Wt and 
5$*(r,t) = 8<S>*{r)e- wt (note that <5$(r) and <5$*(r) are 
independent), and consequently Sn(r,t) = 5h(r)e~ tuJt , 
6rh(r,t) = 8m(r)e~' luJt and 5n/(r, f) = 5nf(r)e~ luJt , one 
finds 

cj<5$(r) = £<J$(r) + g&f) n (r)<5$*(r) 

+2p 6fe $ (r)5n(r) + gbb $ (r)Jm(r) 
+ 56/ $o(r)5n/(r). (7) 

In the absence of coupling terms, Eq. Q and its ad- 
joint are formally equivalent to the time-independent 
Bogoliubov-deGennes (BdG) equations [l3| 



-Wi(r) 



C gbbn (r) 
gb b n (r) C 



Ui(r) 
Vi(r) 



. (8) 

which define the Bogoliubov quasiparticle wave functions 
Mj(r) and Vi(r) with excitation energies ujf . This equiv- 
alence is not surprising since the Bose broken symmetry 
leads quite generally to the one-one correspondence be- 
tween the small oscillations of the condensate and the 
single-quasiparticle wave functions |26| . For the purpose 
of solving Eq. Q , to leading order of the two coupling- 
constants [23, we thus can select the Bogoliubov quasi- 
particle wave functions corresponding to the low-energy 
collective mode that we are interested, and set accord- 
ingly <5$ (r) = u(r), 8$%(r) = v(r) and uj = uj . 

The first-order correction due to the fluctuations 5n(r), 
5m(r) (and its complex conjugate), and <5n/(r), can be 
calculated by expanding 



5$(r) 
<5$*(r) 



uj 



u(r) 
v(r) 



luq + Suj — 17, 



<S$i(r) 



(9) 



(10) 



where Soj represents the shift in the real part of the fre- 
quency and 7 is the damping rate. The correction of 
the wave functions in Eq. J§J) is chosen to be orthog- 
onal to the unperturbed Bogoliubov quasiparticle wave 
functions, 



dr (u* (r) 5$ x (r) - v*{r)5$l{rj) = 0. (11) 



Inserting this perturbation ansatz into Eq. Q and its 
adjoint, wc multiply the first equation by u* (r) and the 
latter by v* (r), and integrate over space. By using Eq. 
(|ll[l and the normalization condition 



dr (u* (r) u(r) — v*(r)v(r)) = 1, 



(12) 



we get the following relation for the eigenfrequency cor- 
rection: 



Suj — 17 



dr<S> (r) [2g b b (u*(r)+v*(r))6n(r) 

+g b b (u*(r)5rh(r) + v*(r)5m*(r)) 

+ g bf (u*(r)+v*(r))6n f (r)}. (13) 

In the next subsection, based on the RPA and linear re- 
sponse theory, we will derive the explicit expressions for 
8h(r), 8m(r), 5m*(r) and <5n/(r), which are induced by 
the condensate oscillations. An alternative way to get 
these expressions in case of pure Bose gases has been 
outlined by Giorgini in Ref. |16| |. Our derivation pre- 
sented below is somewhat simpler and more transparent 
in physics. 



B. RPA and linear response theory 

Let us consider the interaction terms in the density 
Hamiltonian (JIJ that couples the condensate wave func- 
tion to the noncondensate component and the degenerate 
Fermi gas 



Hi 



= ^ / dr 



4|$(r,t)| 2 V + M)V(r,t) 



+$* 2 (r, t)$(r, t)$(r, t) + $ 2 (r, t)$ + (r, t)i> + (r, t) 

+g bf J dr|$(r,i)| 2 0+(r,i)<Kr,i). (14) 

Consistent with setting the triplet averages to zero in 
derivation of Eq. J2J), we have dropped the terms linear 
in $(r, t). In the spirit of RPA, by linearizing the above 
interaction Hamiltonian, we identify the perturbation in- 
duced by small amplitude-oscillations of the condensate 
(the e~ iuJt dependence is not shown explicitly): 

nZt = 9bb |dr$o(r) [2Kr) + ! ;(r))^(r,^(M) 

+v(r)rp(r, t)$(r, t) + u(r)ijj + (r, i)-0 + (r, t) 

+g bf J dr$ (r) (u(r) + v(r)) 4> + (r, t)0(r, 0$ 
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5m 


= 9bb 


\ 5m* j 






where to the leading order we have replaced SQ(r,t) 
and 5$>*(r, t), respectively, by u(r)e~' lujt and v(r)e~ tuJt . 
Within the linear response theory, the fluctuations are 
given by 



Xfirh Xnfh+ 
Xrhrh X r mrh+ 
Xrh+rh Xrh+rh+ 



and 

Sn f = g bf [ dr'xf (r, r'; w) $ (r') (u(r') + v(r')) . (17) 

Here we define Xap^ou = J dr'x a p (r, r'; uS) $ (r')u(r') 
and Xap^ov = J dr' Xal 3 (r,r';w) $ (r>(r') in Eq. JTSJl, 
where the indices a,/3 = n,m, or m + . an d X/ are 
the usual two-particle correlation functions for the Bose 
and Fermi gas 28] . By using Wick's theorem, they can 
be easily expressed in terms of the quasiparticle ener- 
gies and wave functions. For instance, for x™, with 
the help of the Bogoliubov transformation we can write 

■0(r, t) = J2i Ui(r)&ie~ iu * * + v* (r)af e^* ' in terms of 

the Bogoliubov quasiparticle operators dj and af . It is 
then straightforward to obtain 



one finds to second order of g b b and gbf, 
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+ (u*v* + v*u*) (1 + ff + ff) 



or 



-(■ 



k*«* + + (i + /f + ff) 



where u> + = u + i0 + . We have used the abbrevia- 
tion: u^UjitiU^ - ..•^•..f^-.f-'^.*r->-\ fB 



uK^Uji^Uiir^u^r'), etc., and // 
(d^di) = 1/ (eP Ui — ij is the Bosc-Einstein distribu- 
tion function with = l/fc^T. x- 1 ^ an( l XaI correspond, 
respectively, to the excitation of single quasiparticles and 
of pairs of quasiparticles. For x/, we have 



Xf 



y 



'I) 



(19) 



1 ] is the Fermi-Dirac distribu- 



where ff = 1/ (e 

tion, and the single-particle wave function ^j(r) satisfies 
the stationary Schrodinger equation |13| 



V 

2m/ 



+ VtLp 0) - M/ + 56/«b (r) 



F 

Uli (pi. 



(20) 



C. eigenfrequency correction 



Substituting the fluctuations (|16|) and l|17|) into Eq. 
(|13fl . and using the explicit expressions for and Xfi 



E(/f-// 



^+ + (wf - wf ) 



\B, 



\Bi 



ui+ + (wf + cjj 3 ^ 



lev 



' ; w+ + (ujf - wf 



(21) 



where the matrix elements Ay , Bij, B^, and Cy are, 
respectively, given by 



As 



(18) 



B , 



B^ - 



— 



dr$ Q [u (uiUj + v. t v* + ViU*) 

+V [uiU* + ViV* + UiVj)] , 

J dr$ [u (u*v* + v*u* + u*u*) 
+v (u*v* + v*u* + v*v*)] , 
J dr$o [u (uiVj + ViUj + ViVj) 
+v (uiVj + ViUj + UiUj)] , 

dr$ (r) (u (r) + v (r)) ^ (r) ^ (r) . (22) 



Eqs. I|21|) and (|22|l are the main result of this section. 
Without the fermionic component, these equations co- 
incides with the finding obtained by Giorgini (the Eqs. 
(39) and (40) in the second paper of Ref. [l(j ) as they 
should be. The last term in the right-hand side of Eq. 
(|21|l is novel and arises from the many possibilities of 
independent particle-hole excitations [291 ]. This mecha- 
nism is known as Landau damping due to the Bose- Fermi 
coupling. On the other hand, the first and second terms 
in the right-hand side of Eq. I|21() correspond, respec- 
tively, to the Landau and Beliaev processes due to the 
interaction between bosons [T^ . 

One of the advantages of our derivation presented here 
is that to obtain the eigenfrequency correction we don't 
need to impose any constraint used in solving the equi- 
librium problem, i.e., the Popov prescription m (r) = 0. 
For a pure Bose gas at high temperatures close to T c , it 
might be reasonable to use the Hartree-Fock spectrum 
for Xa/3- As a result, only xah is nonzero and the eigen- 
frequency correction reads 



Suj 



i 1 = 4g 2 bb J dr J dv'Qo (r) (u* (r) + v* (r)) 

x X nn (r, r'; u) $ (r') (u (r') + v (r')) (23) 



which is identical to the finding of Reidl et al. obtained 
by using the dielectric formalism (the Eq. (52) in Ref. 
|), if one notices that Sn c (r) = $o ( r ) (u (r) + v(r)). 
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It should be noted that the second term in the right- 
hand side of Eq. (|2*T)l . corresponding to the Beliaev 
process, is ultraviolet divergent. This reflects the fact 
that the contact interaction is an effective low-energy in- 
teraction invalid for high energies. One way to remove 
this divergence is to express the coupling constant gbb in 
terms of the two-body scattering matrix obtained from 
the Lippman-Schwinger equation. This renormalization 
scheme has been put forward in Ref. ,10], however, only 
valid in the thermodynamic limit, where the Thomas- 
Fermi approximation can be implemented [T^ . In our 
derivation, one can explicitly show that such divergence 
comes from the two correlation functions: Xmm+ ( r i r '; w ) 
and Xm+m ( r 7 r 'y UJ )- One thus may wish to remove the 
divergence by regularizing Xmm+ an d Xm+m m real space 
in a way similar to that described in Refs. |3Ct 13 1| . In 
this paper, for simplicity we shall neglect the second term 
in the right-hand side of Eq. H21[l , since it is always very 
small compared with other two terms at all tempera- 
tures. This treatment is well justified by the excellent 
agreement between the experimental result [3^ | and the 
theoretical prediction by Reidl et al. j^j for a pure Bose 
gas, where in the theoretical calculations the Beliaev pro- 
cess is completely ignored. 

There is one last technical issue to resolve: concerning 
the damping rate, the terms in Eq. (|21|) involve a sum 
over many 5 functions in energy, which, if interpreted ex- 
actly, will tend to be null for discrete quasiparticle states. 
Here we shall adopt the strategy of Ref. 33] and use an 
expression with a Lorentz profile factor in place of the 
energy S function 



Suj — i^f = 



-ah X! 



(/f-//)iA,f 



LOQ + Suj + (uf - to?) + i<y 

{f?-ff)\Cg\* 
uj a + Suj + (uf - ujf) + ij 



,(24) 



which can be solved iteratively for Suj and 7. This expres- 
sion can be formally obtained from Eq. 1)21(1 by assuming 
that the perturbed resonance frequency u> is distributed 
over a range of values characterized by a Lorentz profile 
with a width 7. 

The structure of our calculation is then as follows: 
First we solve the unperturbed equilibrium problem for 
ujq , u(r) and v(r). This step requires solving a closed set 
of Eqs. 0, (|SJ), and (|20() . which we have referred to as 
the "HFB-Popov" equations for dilute Bose-Fermi mix- 
tures. We already have reported on our self-consistent 
algorithm in Ref. JlU for this problem. As a result, 
we have all the necessary inputs, namely, $o(r), uf, 
uf , Uj(r) and Wi(r), and <fii(r) for performing the second 
step: the use of Eq. I|24|) in connection with the matrix 
elements Aij and Cy defined in Eq. 



III. NUMERICAL RESULTS 

In this work we analyze the low-lying condensate oscil- 
lations of a Bose-Fermi mixture for varying Bose-Fermi 
coupling constant and temperature in an isotropic har- 
monic trap, for which the order parameter $o( r ), the 
Bogoliubov quasiparticle amplitudes Ui(r) and «i(r), and 
the orbits (pi(v) can be classified according to the num- 
ber of nodes in the radial solution n, the orbital angu- 
lar momentum I, and its projection m. For this sake, 
we apply the theory developed in the proceeding section 
to mixtures of 41 K— 40 K and 87 Rb— 40 K. For the former 
system, we discuss some generic properties of Bose-Fermi 
mixtures. In this case the rigid oscillation of the center of 
mass, or the dipole mode, will also be an eigenstate of the 
many-body system. As a result, the oscillation frequency 
will be fixed at the bare trapping frequency, regardless of 
any interactions. The fulfillment of this property, which 
is usually referred to as the generalized Kohn theorem, 
thus provides us a stringent test on the correctness of our 
results. The second choice of 87 Rb— 40 K mixture corre- 
sponds to a specific example available in present exper- 
iments 0. In this case we include explicitly the mass 
difference and the different oscillator frequencies of the 
trapping potentials for the two species, and build the 
possible experimental relevance of our results. 



A. 41 K- 40 K 

We first consider a mixture of 2 x 10 4 41 K (boson) 
and 2 x 10 4 40 K (fermion) atoms with the following 
set of parameters: nib = mj — 0.649 x 10~ 25 kg , 
ujb = — 2n x 100 Hz, = 286ao = 15.13 nm |34j . 
where ao = 0.529 A is the Bohr radius. We also express 
the lengths and energies in terms of the characteristic os- 
cillator length a b ho — (h/nibUb) 1 ^ 2 and characteristic trap 
energy TiUb, respectively. 

In Figs, (la) and (lb), we present, respectively, our 
results for the monopole (I = 0) and quadrupole (I = 2) 
condensate oscillations at a very low temperature T = 



0.01T C °, where T c ° = 0.94w b W fc 1/3 / k B w 122 nK is the 
critical temperature for an ideal Bose gas in the thermo- 
dynamic limit. The mode frequencies, in units of the bare 
trapping frequency, are plotted as a function of the Bose- 
Fermi coupling constant measured relative to the Bosc- 
Bose coupling constant, k — gbf /gbb- The lines with open 
circles show the unperturbed frequencies obtained from 
the HFB-Popov equations, uio, while the lines with solid 
circles denote the values after correction, uj — ujq + Slu. 
For comparison, the predictions of the scaling theory at 
zero temperature are also plotted by the dashed lines 
[8j- At this low temperature, the eigenfrequency shift 
Suj arises mainly from the dynamics of the degenerate 
Fermi gas. For small values of |k| < 1, Suj is negligibly 
small due to square dependence on the Bose-Fermi cou- 
pling constant. However, as |k| increases Suj becomes re- 
markable. In particular, the corrected frequency for the 
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FIG. 1: The dispersion relation of the monopole (I = 0) and 
quadrupole (I = 2) condensate oscillations for a mixture com- 
posed of 2 x 10 4 41 K (boson) and 2 x 10 4 40 K (fermion) atoms 
at T = O.OIT^ 1 , where T c ° = 0.94w fc AT i 1/3 /k B ~ 122 nK is 
the critical temperature for an ideal Bose gas in the thermo- 
dynamic limit. The mode frequencies, in units of the bare 
trapping frequency, are plotted as a function of the reduced 
Bose- Fermi coupling, k = gtf/gtb- The lines with open cir- 
cles show the unperturbed frequencies calculated by the static 
HFB-Popov equations &o, while the lines with solid circles de- 
note ui = ujo + Suj. For comparison the predictions of the scal- 
ing theory in Ref. [8] are also plotted by the dashed lines. The 
inset in (b) shows the dispersion relation of the dipole excita- 
tion. The other parameters used in the numerical calculation 
are: m b = m f = 0.649 x 1CT 25 kg, uj b = Uf = 2iv x 100 Hz, 
and abb = 286ao = 15.13 nm, where ao = 0.529 A is the Bohr 
radius. 



FIG. 2: The same as in FIG. 1, but for T = 0.75T C °. In 
the insets the solid and dashed lines show, respectively, the 
fractional shift Suj/ljo due to the first and second terms in Eq. 

Eg . 
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quadrupole mode decreases with increasing Bose-Fermi 
coupling constant up to k « 3, at which a sharp upturn 
occurs. This sharp dip is accompanied by a dramatic in- 
crease of damping rates (not shown in the figure), and 
can be well interpreted as a signal to approach the spa- 
tial separation (demixing) point of the two species. On 
the contrary, the unperturbed quadrupole frequency loq 
has a qualitatively different behavior against k: it shows 
a parabolic dependence with a minimum located around 
K = 0.5. In spite of its large value at |«| > 1, Suj is still 
much smaller than unperturbed frequency ujo for all the 
calculated points. Therefore, the criterium for the appli- 
cability of the perturbation theory in Eqs. @ and 1)10(1 is 
justified. In the inset of Fig. (lb) we also show the result 
for the dipole mode (1 = 1). As we can see, although ujq 
deviates significantly from the bare trapping frequency 
u>b for relatively small values of k, uj still persists at ujb 
over a wide range of k (i.e., —2 < n < +3). This is con- 
sistent with the generalized Kohn theorem. As a result, 
the correctness of our theory and numerical calculations 
is partly checked. 

The eigenfrequency shifts are affected by the temper- 
ature. In Fig. 2, we report the mode frequencies against 
k at a high temperature T = 0.75T°, where the conden- 
sate oscillates in the presence of a large fraction of above- 
condensate atoms. Compared with the results for the low 
temperature, the eigenfrequency shifts are considerably 



FIG. 3: The damping rates of the monopole (a) and 
quadrupole condensate oscillations (b) at T = 0.75T,?. The 
solid and dashed lines correspond to the contribution from 
the Landau process due to the Bose-Bose interaction and due 
to the Bose-Fermi coupling, respectively. The line with solid 
circles is the total contribution. 

reduced. The sharp dip at n « 3 for the quadrupole mode 
also becomes much broader. Moreover, in the absence 
of the Bose-Fermi coupling, the eigenfrequency shift is 
nonzero for the quadrupole mode. This is caused by the 
dynamics of the noncondensate component as shown in 
the inset of Fig. (2b), where the solid and dashed lines 
depict, respectively, the fractional shift Soj/ujq due to the 
Bose-Bose interaction and Bose-Fermi coupling (or, in 
other words, due to the first and second terms in Eq. 

The damping rate of condensate oscillations at finite 
temperature deserves its own study. In Figs. (3a) and 
(3b), we show, respectively, our predictions on the damp- 
ing rates of monopole and quadrupole oscillations at 
T = 0.75X 1 ,? . The lines with solid circles are the sum 
over two contributions: one is the Landau damping due 
to the Bose-Bose interaction (the solid lines), jbbi an d 
the other is the Landau damping due to the Bose-Fermi 
coupling (the dashed line), jbf- For the monopole mode, 
the essential feature is the decrease of the damping rate 
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FIG. 4: The damping strength 7^ defined in Eq. 1261 as 
a function of the transition frequencies (in units of u>b) 
allowed by the selection rules, (a) 7^ for the monopole ex- 
citation and (b) j^j for the quadrupole excitation. In each 
subplot, upper and lower panels correspond to the case of 
k = —2 and k = +6, respectively. In the former case the 
mixture is in mixed regime, while in the latter the bosonic 
and fermionic density profiles separate in space. The arrows 
point to uj — wo + Suj. 



as k increases towards the demixing point. This decrease 
is mainly attributed by jtb, and reflects the reconstruc- 
tion of the bosonic monopole-excitation spectrum across 
the demixing point. To better illustrate this point, we 
rewrite the expression for 7 in the following form 



7 
-fbb 

Ibf 



Jbb + Jbf, 



(25) 



7 » 



7/tt 



w + Suj + (wf - ujf)] + 7 2 



bf 



7/tt 



[u) + 5uj+ (uj[ - uf)] " + 7 2 



where the "damping strength" 



l bb = 



-|^| 2 (/f-/f) 
UjQ 



(26) 



have the dimensions of a frequency. In Fig. (4a), we plot 
7y against the transition frequency = ujf — ujf > 
allowed by the selection rules for n = —2 and k = 6. For 
the latter value of k, the overlap between the bosonic 
and fermionic cloud is very small, and the mixture is 
deep into the demixing regime. Compared with the mix- 
ing case of k = —2, the region of transition frequencies 
at k = 6 narrows, and its center moves to the low-energy 
side. Contrarily the calculated u>o + 5ui has a blue shift 
and is completely out of the transition region. As a con- 
sequence, the condensate oscillation is not damped by the 
Landau process due to the Bose-Bose interaction. For the 
quadrupole mode, instead we observe that the damping 
rate increases as the mixture moves towards the demix- 
ing point with increasing k. This trend comes from the 
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FIG. 5: The dispersion relation of the monopole (a) and 
quadrupole excitations (b) against the reduce temperature 
T/T° at k = —2. The inset in (a) shows the dispersion rela- 
tion for the dipole mode. 




TIT 



TIT 



- c c 
FIG. 6: The same as in FIG. 5, but for the damping rates. 



increase of jbf, and reflects, on the other hand, the re- 
construction of the fermionic quadrupole-excitation spec- 
trum. As shown in Fig. (4b), with increasing Bose- Fermi 
coupling the damping strength 7^/ becomes larger and 
denser. Accordingly, the condensate oscillation is heav- 
ily damped by generating many particle- hole excitations. 

The last study in this subsection concerns the temper- 
ature dependence of the eigenfrequency shifts and damp- 
ing rates at a specific Bose- Fermi coupling constant. In 
Figs. (5a) and (5b), we report, respectively, our results 
for the monopole and quadrupole mode frequencies as 
a function of the reduced temperature T/T° at n = —2. 
The corresponding damping rates are shown in Figs. (6a) 
and (6b). The fulfillment of the generalized Kohn theo- 
rem is checked in the inset of Fig. (5a), where the cal- 
culated dipole frequency is very close to uib (or, more 
precisely, 0.97 < Lo D jui b < 1.0) for all the temperatures 
considered. At this specific value of k, one can see that 
both the monopole and quadrupole frequencies have a 
downshift with increasing temperature, analogous to the 
results obtained for a pure Bose gas jla . 120] . In addition, 
the behavior of the damping rates is also qualitatively 
similar USUI Hi. 
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FIG. 7: The dispersion relation of the monopole and 
quadrupole condensate excitations for a mixture consisting 
of 2 x 10 4 87 Rb and 2 x 10 4 40 K atoms at T = 0.75T C °, where 
T° ~ 112 nK. The other parameters used in the calculation 
are: m b = 1.45xl0" 25 kg, cu b = 2ttx91.7Hz, m f /m b = 0.463, 
uif/ub — 1.47, and abb = 99ao = 5.24 nm. 
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FIG. 8: The same as in FIG. 7, but for the damping rates. 
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We now turn to consider a 87 Rb— 40 K mixture com- 
posed of 2 x 10 4 bosonic and 2 x 10 4 fermionic atoms 
under the conditions appropriate to the LENS experi- 
ments As in experiment, we introduce the quantities 
a = mf/mb = 0.463 and (3 = LUf/ub = 1-47 to parame- 
terize the different mass and different trapping frequency 
of the two species, which satisfy the constraint a(3 2 = 1 
since both bosons and fermions experience the same trap- 
ping potential. In addition, we take the s-wave Bose-Bose 
scattering length abb = 99ao = 5.24 nm j3||, and fix the 
trapping frequency tub = 2tt x 91.7 Hz, which is the geo- 
metric average of the axial and radial frequencies of Ref. 
0. The s-wave Bose- Fermi scattering length is varying, 
and in the experiment it can be conveniently tuned by 
the Feshbach resonance j^. Notice that the calcula- 
tions presented here are restricted to the isotropic traps, 
opposite to the cylindrical symmetric traps used in ex- 
periments. As a result, our results are only useful in a 
qualitatively level. 

In Fig. 7, we plot the frequencies for the monopole and 
quadrupole oscillations as a function of k at T — 0.75X 1 ,?, 



where T 1 ? « 112 nK. Both the monopole and quadrupole 
frequencies decrease slowly with increasing k up to k w 5. 
Above this value the frequencies gradually rise up. In 
the whole region of k, the variation of frequencies due to 
Bose-Fermi coupling is small. However, it is still possible 
to be detected by the accurate frequency measurement. 
For instance, at k — —6, we find that the values of the rel- 
ative variation (ui K= -e — ^k=o) /^k=o for the monopole 
and quadrupole modes are, respectively, 2.8% and 7.6%, 
well within the experimental resolution. 

The damping rate of the monopole and quadrupole 
modes at the same temperature is shown in Figs (8a) 
and (8b), respectively. The behavior of the damping rate 
against k, that is, the decrease (increase) of the monopole 
(quadrupole) damping rate across the demixing point 
k « 5, is very similar to that in Figs. (3a) and (3b), 
except that the overall magnitude is two times smaller. 
This behavior together with the slow rise up of the mode 
frequency around k « 5 thus may provide us a useful sig- 
nal to locate the onset of the phase transition of spatial 
separation. 



IV. CONCLUDING REMARKS 

In this paper we have developed a theory for study- 
ing the low-lying condensate oscillations of a spher- 
ically trapped Bose-Fermi mixture at finite tempera- 
ture in the collisionless regime. In this theory, the un- 
perturbed mode frequency is firstly calculated within 
the static Hartree-Fock-Bogoliubov-Popov approxima- 
tion. The frequency correction, arising from the coupled 
dynamics of the condensate, noncondensate, and degen- 
erate Fermi gas, is then taken into account perturba- 
tively by means of the random phase approximation. We 
have applied our theory to the mixtures of 41 K— 40 K and 
87 Rb— 40 K, and have studied the dispersion relation of 
the monopole and quadrupole condensate excitations as 
a function of the Bose-Fermi coupling at various tem- 
peratures. The correctness of our theory and numeri- 
cal calculations is partly checked by the fulfillment of 
the generalized Kohn theorem for the dipole excitation. 
At a relatively high temperature we find that, as the 
mixture moves towards demixing point with increasing 
Bose-Fermi coupling, the damping rate of the monopole 
(quadrupole) excitation increases (decreases). This be- 
havior provides us a possible signature to identify the 
phase transition of spatial separation. 
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